PSRCHIVE and psrfits: Definition of the Stokes Parameters and 
Instrumental Basis Conventions 

W. van Straten A ' c , R. Manchester® , S. Johnston 3 , and J. Reynolds® 

A Swinburne University of Technology, PO Box 218, Hawthorn VIC 3122, Australia 
B CSIRO Astronomy and Space Science, PO Box 76, Epping NSW 1710, Australia 
c Email: vanstraten.willem@gmail.com 



Abstract: 

This paper defines the mathematical convention adopted to describe an electromagnetic wave and its 
polarisation state, as implemented in the PSRCHIVE software and represented in the psrfits definition. 
Contrast is made between the convention that has been widely accepted by pulsar astronomers and the 
IAU/IEEE definitions of the Stokes parameters. The former is adopted as the PSR/IEEE convention, 
and a set of useful parameters are presented for describing the differences between the PSR/IEEE 
standard and the conventions (either implicit or explicit) that form part of the design of observatory 
instrumentation. To aid in the empirical determination of instrumental convention parameters, well- 
calibrated average polarisation profiles of PSR J0304+1932 and PSR J0742— 2822 are presented at 
radio wavelengths of approximately 10, 20, and 40 cm. 
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1 Introduction 

Since the early days of pulsar astronomy, observations 
of the polarisation properties of their radio signal have 
yielded a diverse range of information about these stars 
and their environment, includin g the properties of the 
pulsar emission mechanism (e.g. Taylor et al.lll97ll : lLvne fc 



l!988l : lEdwards fc Stappersjl2004l), the geometry of t he 



pulsar magnetospher e fe.g. lRadhakrishnan fc Cookdll969l : 
I Johnston et aLlfeoOjjl ). a nd the properties o f the Galac- 
tic m agnetic field (e.g. Manchester] Il974l : Irian et ail 
2006). Polarimetry also provides additional informa- 
tion th at may be exploit ed for high-precision pulsar 
timing l|van Stratenll200rj ). 

Pulsars are typically observed either at the focus of 
a single antenna or at the phase centre of an antenna 
array. In either case, there are effectively two receptors 
that ideally respond maximally to orthogonal senses of 
polarisation. The statistical properties of the voltages 
induced in the receptors are computed as a function of 
topocentric pulse phase, producing a quantity known 
as the polarisation profile. Well-calibrated polarisation 
profiles require observations of one or more calibrator 
signals and a mathe matical description of the instru- 
mental response (e.g. Ivan Stratenl l2004), 

Even after calibration, there is often confusion re- 
garding the definition of left-handed and right-handed 
circular polarisation (LCP and RCP) as well as the 
sign of Stokes V and/or the position angle of the linear 
polarisation. The uncertainty arises from the different 
conventions that have persisted in the classical physics 
and engineering literatu re. The earliest p ulsar polari- 
sation observations (e.g. Manchester 1971) established 
a convention that is consistent with the Institute of 
Electrical and Electronics Engineers (IEEE) definition 
of LCP and RCP and the definition of Stokes V by 
iKrausi (| 19661 ). This convention differs from the one 



later adopt ed by the International Astronomical Union 
(|IAUlll974l 1. 

The PSR/IEEE convention is implemented by the 
PSRCHIVE softwarcQ and encoded in the psrfits file 
format, both of which have been openly developed in 
an effort to facilitate the exchange of pulsar astronom- 
MsaActiatecibctween observatories and research groups 
(jHotan et al.l [2004). The PSRCHIVE software imple- 
ment s innovative meth o ds of polarisation calibration 
(e.g. Ivan Stratenl |2004 120061) and Faraday rotation 
meas ure estimation fe.g. lHan et al.ll2006l : lNoutsos et al.l 
2008). To properly utilise these algorithms, data must 
be input in a form that complies with the PSR/IEEE 
convention. 

In the following section, the PSR/IEEE convention 
is defined and contrasted with the IAU/IEEE defini- 
tions of the Stokes parameters. Section 3 introduces 
a comprehensive set of parameters that can be used 
to document the differences between an instrumental 
design and the PSR/IEEE convention. In Section 4, 
polarisation profiles of PSR J0304+1932 (B0301+19) 
and PSR J0742-2822 (B0740-28) are provided as a 
point of reference for astronomers aiming to conform 
to the PSR/IEEE convention. 

2 Definitions 

Consider a quasi-monochromatic electromagnetic wave 
with mean frequency u), represented at the origin by 
the transverse electric field vector 



e(t) 



eo 
ei 



ao(t)expi[cj>o(t) + wt] 
ai (t) exp i[4>i(t) + u)t] 



(1) 



Note that the complex argument increases linearly 
with time; this sign convention is commonly encoun- 
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tered in engineerin g texts re.g. |Papoulislll962| : |BraceweUl 
Il965l ; iKrausI Il966l ) and is implicit in the definition of 
most forward discr ete Fourier transform (DFT) im- 
plementations (e.g. Frigo & Johnson 2005). It is also 
adopted in a semin al series of papers on r adio polari- 
metr i c calibration llHam aker et al.l 1 19961 : ISault et al.l 
1996; Hamakc r fc Bregmanl 11999 ). Given the above 
definition, time delays correspond to negative values 
of the phase (j>. 

The polarisation of an electromagnetic wave is de- 
scribed by the second-order statistics of e, as rep- 
resented using the complex 2x2 coherency matrix 
(|Born fc WolJll980l ) 

p= ( e®eW ; e °i; ; eo i! V ( 2) 

N ' \ (eie ) (eiei) J 

Here, the angular brackets denote an ensemble aver- 
age, ® is the matrix direct product, and is the Her- 
mitian transpose of e. The coherency matrix may be 
expressed as a linear combination of Hermitian basis 
matrices, 

3 

p = \ Yl, Sk(Tk = (s ° 0-0 + s ■ °")/ 2 ' (3) 

k=0 

where So is the total intensity, S — (Si, S%, Ss) is the 
Stokes polarisation vector, ctq is the 2x2 identity ma- 
trix, and cr = (cri, CT2, CT3) is a three-vector with com- 
ponents equal to the Pauli matrices. The four Hermi- 
tian basis matrices, consisting of the identity matrix 
and the three Pauli matrices, are defined as 

a ° = ( 1 ) ffl = ( -1 ) 

The Stokes parameters are the projections of the co- 
herency matrix onto the basis matrices, given by 

S k =tT(a k p), (5) 

where tr() is the matrix trace operator. 

Linear transformations of the electric field are rep- 
resented by complex 2x2 Jones matrices. Under the 
operation, e' = Je, the coherency matrix is subjected 
to a congruence transformati on, p' = JpJ . Using the 
axis-angle parameterisation (Britton 2000), t he po lar 
decomposition of a Jones matrix lllamnkorl l200Ch is 
expressed as 

J = J B^OS)^), (6) 

where J = |J| 1//2 , jJj is the determinant of J, B T f l (/3) 
is positive-definite Hermitian, and Rft(<£) is unitary; 
both B T f l (/3) and R^(cj>) are unimodular. Under the 
congruence transformation of the coherency matrix, 
the Hermitian matrix 

Bfh(P) = a o cosh j3 + rh • cr sinh j3 (7) 

effects a Lorentz boost of the Stokes four-vector along 
the rh axis by a hyperbolic angle 2/3. As the Lorentz 



transformation of a spacetime event mixes temporal 
and spatial dimensions, the polarimetric boost mixes 
total and polarised intensities, thereby altering the de- 
gree of polarisation. In contrast, the unitary matrix 

Rri (<W = ""o cos cj> + in ■ cr sin cf> (8) 

rotates the Stokes polarisation vector about the n axis 
by an angle 2(f). As the orthogonal transformation of a 
vector in Euclidean space preserves its length, the po- 
larimetric rotation leaves the total intensity and degree 
of polarisation unchanged. 

2.1 The IAU/IEEE Convention 

As summarised by Hamaker & Bregman (1996) and 
depicted in Figure 1, the IAU/IEEE definitions of the 
Stokes parameters are based on a right-handed Carte- 
sian coordinate system, in which the plane wave prop- 
agates toward the observer in the positive z direction, 
and eo = e x and ei = e y are the components of the 
electric field projected onto North and East, respec- 
tively. Noting that the sign of the complex argument in 
equation (JTJ corresponds to the upper sign convention 
i n equ ations (2) through (4) of Ham aker fc Bregmanl 
(1996), the Stokes parameters are defined by 



I = (\e x \ 2 + \e y \ 2 ) (9) 

Q = (\e x \ 2 -\e y \ 2 ) (10) 

U = (2Re[e x e;]) (11) 

V = <2Im[e x e;]} (12) 



Stokes V is positive for RCP, for which e y lags e x and, 
looking toward the source, the observer sees the elec- 
tric field rotate counter-clockwise. The position angle 
of the linearly polarised component is also measured in 
a counter-clockwise direction from North toward East. 

2.2 The PSR/IEEE Convention 

In the majority of published pulsar polarisation ob- 
servations, Stokes V is positive for LCP as defined 
by the IEEE fe.g | Manchesterl Il97ll : IStinebring et all 
1 19841 : IXilouris et al.lll998l ). This is the convention de- 
scribed bv lKrausl (|l966T ): it is encoded in the sign of the 
phase in equation ([l]), the definition of the Pauli matri- 
ces in equation @, and the relationships between the 
Stokes parameters and the coherency matrix defined 
by equations ((3} and (5). Using either equation © or 
equation © with S = (Q, U, V), it can be shown that 
in the linear basis, 

n 1 ( I + Q U-iV \ ( . 

P =2{ U + iV I-Q )■ ^ 

In this case, 

U + iV = 2{e y e* x ). (14) 

That is, Stokes V is positive when the phase of e y 
leads that of e x , which is opposite to the IAU con- 
vention for Stokes V. The PSR/IEEE convention is 
adopted and used in the remainder of this paper; it is 
the standard implemented by the psrchive software 
and expressed in the psrfits definition. 
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Figure 1: A right-hand circularly polarised (RCP) wave in the adopted coordinate system. The x and y 
axes lie in the plane of the sky, with the x-axis pointing North, the y-axis pointing East, and the z-axis 
pointing toward the observer. The electric field vector e lies in the x-y plane; its position angle \& is 
measured with respect to North, increasing toward East. As the RCP wave propagates in the positive z 
direction, the position angle increases; that is, the electric field vector rotates counter-clockwise, as seen 
by the observer. Note that an RCP wave, as defined by the IEEE, forms a left-handed helix in space. 



The PSR/IEEE convention for the position an- 
gle of the linear polarisation is consistent with the 
IAU/I EEE convention. However , as discussed in de- 
tail bv lEverett fc Weisberd (120011). the rotatin g vector 
model fRVM; iRadhakrishnan fc Cookd Il969f ) has of- 
ten been applied to measurements using an equation 
in which the position angle is assumed to increase in a 
clockwise direction on the sky, which is opposite to the 
IAU/IEEE convention. Therefore, care must be taken 
when comparing RVM parameter estimates from dif- 
ferent experiments. 

Given the definition of the Stokes parameters in 
terms of linearly polarised receptors, it is possible to 
also define a standard basis for circularly polarised 
feeds. The basis transformation, C, is chosen to effect 
a convenient cyclic permutation of the Stokes param- 
eters; that is, under the congruence transformation, 
p' = CpC\ the Stokes vector, S' = (V, Q, U). 

The transformation with unit determinant that ef- 
fects the desired cyclic permutation of the Stokes pa- 
rameters is 



1 -i 
1 i 



= {e L e R )\ 



(15) 



where 



and en 



1 

—i 



are the Jones vectors of left and right circularly po- 
larised receptors expressed in the linear basis. Apart 
from an absolute phase factor of e"" 1 ^ 4 chosen to yield 
a basis transformation with unit determinant, the above 
definition of circularly polarised recep tors is identical 
to that of Ha maker k, Bregmanl (|l996h . 



3 Instrumental Conventions 

A wide variety of different conventions, both explicit 
definitions and implicit assumptions, are encountered 
in the design of observatory instrumentation and data 
analysis software. Therefore, it is desirable to define a 
set of parameters that describe any known differences 
between the experimental design and the chosen stan- 
dard. These corrections can be be broadly classified 
into four groups: projection, basis, instrument, and 
backend. 

The known projection of the ideal feed receptors 
onto the plane of the sky may include effects such as 
the parallactic rotation of an altitude-azimuth mounted 
antenna, or the foreshortening of fixed dipoles in an ar- 
ray. The ideal basis describes a feed consisting of two 
orthogonally polarised receptors with unit gain. The 
feed may also be equipped with an internal reference 
source, such as a noise diode, that may be used to 
calibrate the differential gain and phase of the instru- 
mentation. The instrument parameterises the trans- 
formation from the receptors to the backend, as de- 
termined using a separate calibration procedure (e.g. 
van Stratenl 120041 ). The backend includes all hard- 
ware and software involved in reducing the astronomi- 
cal data to produce results that are suitable for further 
analysis. Only the basis and backend correction pa- 
rameterisations are discussed here. 



3.1 Basis Parameters 

In principle, the polarisations of both an ideal feed 
and the reference source could be specified by four 
rotatio ns; e.g. two sets of el lipticity and orientation 
angles ( Chandrasekhar 1960). However, basis rota- 
tions can be problematic due to singularities at the 
poles. Furthermore, it may be preferable to use less 
abstract, more commonly encountered, and/or more 
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readily measured quantities. The convention may also 
be simplified by assuming that the receptors are either 
linearly or circularly polarised, and that the reference 
source is linearly polarised. Under these assumptions, 
the following four parameters are sufficient to com- 
pletely describe the receiver. 

Feed Basis: circularly or linearly polarised receptors 

Feed Hand: left or right; if left, eo and ei are switched 

Symmetry Angle: orientation of a linearly polarised 
signal that produces an equal and in-phase response in 
each receptor 

Calibrator Phase: differential phase of the internal 
reference source, $ c = tan^ 1 (S3 / S2) 

In the linear basis, a linearly polarised wave os- 
cillating in the Northeast-to-Southwest plane (positive 
Stokes U) will produce equal and in-phase responses 
in each receptor; therefore, the symmetry angle has a 
nominal value of 45 degrees. Also, because the refer- 
ence source is linearly polarised (S3 = 1^ = 0), there 
are only two possible values for the calibrator phase: 
(U > 0) and 180 (U < 0) degrees. 

In the circular basis, a linearly polarised wave os- 
cillating in the North-to-South plane (positive Stokes 
Q) will produce equal and in-phase responses in each 
receptor; therefore, the symmetry angle has a nominal 
value of degrees. Also, because S2 = Q and S3 = U, 
the calibrator phase is arbitrary. 

Given these parameters, the Jones matrix of the 
feed is given by 

J = XCR C (6), (16) 

where X is the identity in a right-handed system and 
the exchange matrix in a left-handed system, C is the 
identity in the linear basis and given by equation (|15[1 
in the circular basis, and B is equal to the symmetry 
angle less its nominal value (45 degrees in the linear 
basis and degrees in the circular basis). Note that, 
as defined by equation P|. R#(0) effects a rotation of 
the Stokes polarisation vector about the Stokes V axis 
by an angle 20. The coherency matrix of the internal 
reference source is 

_ 1 ( 1 exp(-z<I> c ) \ 

Pc ~ 2 ^ exp(i$ c ) 1 ) ' [U) 

where $ c is the calibrator phase. 

3.2 Backend parameters 

To describe the design of the backend, it is sufficient to 
consider only the complex conjugation of the electric 
field. This occurs during lower sideband down con- 
version and/or when the design of an instrument is 
based upon a different convention for the sign of the 
phase than that adopted in equation fTJ. in the case of 
dual-sideband down conversion, complex conjugation 
occurs when the phase of the quadrature component 
leads that of the in-phase component of the signal. 
Complex conjugation causes a sign change in S3, a re- 
flection that cannot be represented by a Jones matrix. 



The following two parameters are used to describe the 
backend. 

Backend Phase: positive or negative 

Down Conversion: true or false; if true, the con- 
jugation due to lower sideband down conversion has 
been corrected 

The correction of complex conjugation is performed 
when either of the following conditions is exclusively 
true: 

A) backend phase is negative, or 

B) down conversion is false and lower sideband down 
conversion was used. 

The correction is not performed when both A and B 
are true because they negate each other. 

3.3 Implementation 

Each of the parameters described in the preceding sec- 
tions has a corresponding parameter in the psrfits 
definition and a representation in the psrchive soft- 
ware. Table Q] summarises these parameters, the range 
of acceptable values for each parameter, and the effects 
that they have on calibrated data. 

4 Discussion 

In general, it is non-trivial to determine the instru- 
mental convention parameters from first principles and 
an empirical approach must be adopted. This typi- 
cally requires observations of either artificial sources 
with known polarisation or astronomical sources that 
have been previously well calibrated, such as the aver- 
age polarisation profiles of PSR J0304+1932 and PSR 
J0742— 2822 presented in Figure 2. These sources are 
visible to observatories in both Northern and Southern 
hemispheres; they also have relatively high flux den- 
sities, well-defined position angle variation, and suffi- 
cient circular polarisation for the purposes of compar- 
ison with other observations. 

These data were recorded at the Parkes Observa- 
tory on 2009 October 9 and 10 using the 64 m antenna 
and two of the Parkes Digital Filter Banks (PDFB3 
and PDFB4). For each pulsar, a 256 MHz band cen- 
tred at 1369 MHz was observed using the centre beam 
of the 20 cm Multibeam receiver and two bands, one 
with 64 MHz of bandwidth centred at 732 MHz and 
one with 1024 MHz of bandwidth centred at 3.1 GHz, 
were observed using the 10-50 cm dual-band receiver. 
Prior to each pulsar observation, a noise diode coupled 
to the receptors was observed and used to calibrate 
the differential gain and phase of the instrument. For 
the 20 cm observations, the receptor cross- coupling pa- 
ramet ers were calibrated as described in Ivan Stratenl 
(|2004T ). Absolute flux calibration was performed using 
observations of the noise diode recorded while point- 
ing on and off the bright Fanaroff-Riley type I radio 
galaxy, 3C 218 (Hydra A), which was assumed to have 
a flux density of 43.1 Jy at 1400 MHz and a spectral 
index of 0.91. 
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Figure 2: Average polarisation profiles of PSR J0304+1932 (left) and PSR J0742-2822 (right), observed 
at radio wavelengths of approximately 40cm (top), 20cm (middle), and 10cm (bottom). In each plot, 
the lower panel displays the total intensity, Stokes / (solid); the linearly polarised flux, L = \J Q" 1 + U 2 
(dashed); and the circular polarisation, Stokes V (dotted); the upper panel plots the position angle of the 
linear polarisation, W — 0.5 tan (U/Q) (the error bars are two standard deviations in length). Stokes V 
is positive for left-hand circularly polarised radiation as defined by the IEEE. 
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Table 1: Instrumental convention parameters and their effects in each basis. 



Parameter 


PSRFITS 


Range 


Effect 




Name 




Linear 


Circular 


Backend Phase 


BE.PHASE 


±1 


± V 


± U 


Down conversion 


BE.DCC 


0/1 


± V 


± u 


Feed Basis 


FD.POLN 


LIN or CIRC 


(Q,u,v) 


(V,Q,U) 


Feed Hand 


FD.HAND 


± 1 


± Q&V 


± U&V 


Symmetry Angle 


FD.SANG 


-tt/2 < 6 < tt/2 


R*(0-7r/4) 


R*(0) 


Calibrator Phase 6 


FD.XYPH 


— 7T < $ c < 7T 


± u&v 


R«($ c /2) 



Applies only when down conversion is effectively lower sideband. 
b In the linear basis, <I> C = or ir. 



The data were corrected for Faraday rotation us- 
ing the most recently published rotation measure (RM) 
for each pu lsar, as determined usin g the ATNF Pulsar 
Catalogue dManchester et alJ2005f). F or PSR J0304+1932, 
RM=-8.3rad m~ 2 i|Manchesterl 1971); for P SR J0742-2822 
RM=+149.95radm 2 l|Johnston et al.ll2005h . The po- 
sition angles plotted in Figure 2 are referred to the 
centre frequency of the band and are as observed at 
Parkes. Variations in the Faraday rotation along the 
line of sight, particularly those occurring in the Earth's 
ionosphere, will rotate the observed angles by a small 
amount. The ionospheric RM component was esti- 
mated using the Int ernational Reference Ionosphere 
models (|Bilitzall200iT ) and was about —0.7 rad m 2 in 
the direction of PSR J0304+1932 and -1.5 rad m" 2 
in the direction of PSR J0742-2822. Position angles 
above the ionosphere were therefore greater than the 
plotted values by about 0?4 and 0?8 at 3100MHz, 0?7 
and 1?5 at 1369 MHz, and 7? 5 and 14? 4 at 732 MHz for 
PSR J0304+1932 and PSR J0742-2822 respectively. 

Referring to the effects defined in Table 1, the 
instrumental convention parameters of new systems 
can be adjusted until the observed values of Stokes 
V and position angle match those plotted in Figure 2. 
This comparison should be performed only after the 
projection (e.g. parallactic rotation) and instrument 
(e.g. differential gain and phase) have been calibrated. 
Given the overlapping domains of the parameter ef- 
fects, the proper selection of parameter values will 
most likely require consultation with both engineering 
documentation and observatory staff. 
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